Noise correlations in the expansion of an interacting ID Bose gas from a regular array 
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We consider the one dimensional expansion of a system of interacting bosons, starting from a 
regular array. Without interactions the familiar Hanbury Brown and Twiss effect for bosons gives 
rise to a series of peaks in the density-density correlations of the expanded system. Infinitely 
repulsive particles likewise give a series of dips, a signature of the underlying description in terms of 
free fermions. In the intermediate case of finite interaction the noise correlations consist of a set of 
Fano resonance lineshapes, with an asymmetry parameter determined by the scattering phase shift 
of a pair of particles, and a width depending on the initial momentum spread of the particles. 



The Hanbury Brown and Twiss (HBT) effect pQ is a 
fundamental signature of quantum statistics appearing in 
quantum optics, atomic and mesoscopic physics, and nu- 
clear collisions [2HB]. It is most dramatically manifested 
as an interference effect in the intensity correlations due 
to two or more incoherent sources, with a sign depend- 
ing on the statistics of particles: positive correlations for 
bosons; negative for fermions. 

In most known instances of the HBT effect interac- 
tions between particles do not play a significant role, ei- 
ther because these effects are weak or due to the spatial 
separation of the sources. In this Letter we consider the 
one-dimensional expansion of a system of particles, where 
strong interaction effects are unavoidable. Indeed, in ID 
the trajectories giving rise to the HBT effect must cross. 

The situation that we will consider is illustrated in 
Fig. [I] Particles are initially confined to a regular ID 
lattice of spacing A, with one particle per site. At time 
t = the lattice potential is removed, though the po- 
tential restricting the particles' motion to one dimension 
remains. We are concerned with the density correlations 
present after some time t, when the system has expanded 
to many times its original size (analogous to the 'far field' 
limit in optics). Thus we have in mind a ID version of the 
experiment of Ref. [7 , in which noise correlations were 
measured in the expansion of a 3D atomic Mott insulat- 
ing state from an optical lattice. A recent experiment 
demonstrated the preparation of such a ID state in a 
slightly different context [8]. 

To introduce some ideas and notation we briefly de- 
scribe the familiar HBT effect in this setting. We as- 
sume Gaussian initial wavefunctions corresponding to 
harmonic oscillator length t = ^h/muj, (f a (y) = 

( 7r ^2 1 )i/4 ex P — ^ V ~2t^ • The ovei> l a P e -A2 / 4 ^ 2 between 
neighboring sites is assumed to be negligible. After a 
period t of free evolution these wavefunctions have the 
form 
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FIG. 1. ID expansion of atoms from an optical lattice. Noise 
correlations will be present in an absorption image of the ex- 
panded cloud. 



sider a pair of identical particles on sites a and a + 
1, the two-particle wavefunction is ^2(^1,^2;^) = 
-^[(p a (x 1 ;t)(p a+1 (x 2 ;t) ±(p a (x 2 ;t)(p a+1 (x 1 ;t)}, with ± 
for bosons and fermions respectively. The corresponding 
probability density is then 
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where the variables ^i j2 = #1,2 A correspond to the ve- 
locities of the two particles. The oscillatory second term 
describes the HBT effect, with a sign dependent on the 
statistics of the particles. For an array of N particles 
the density-density correlation function develops peaks 
due to the contributions of higher harmonics arising from 
pairs of particles separated by multiples of A 

C(x 1 ,x 2 ]t) = J dx 3 • • 'dx N \V N (x 1 ,x 2 , • • • ,x N ]t)\ 2 . 
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In a trajectory picture the HBT effect arises as a cross- 
term between trajectories that do and do not exchange 
pairs of particles (see Fig. |3j bottom) 
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FIG. 2. (Top) Normalized correlation function 

N ( C ^;g. 2 t f - l) for cA = 2, £/A = 0.2. (Bottom) A 
slice with x\ — —X2 for the same parameters, showing the 
evolution of the Fano asymmetry between successive peaks. 

We turn now to the central subject of this paper: the 
HBT effect in the presence of interactions between the 
particles. We assume that the evolution of the system 
for t > is governed by the iV-particle Hamiltonian 

1 N d 2 

i=l 1 i<j I 



The c — >• and c — >• oo limits can be described in terms of 
free bosons and free fermions, respectively. The density- 
density correlations reflect this, corresponding to the plus 
sign in Eq. ([3| in the former case and the minus sign in 
the latter. Our main result, valid when e~ 2cA <C 1, is 
that in the crossover regime the density-density correla- 
tions consist not of a series of symmetric peaks or dips 
but rather of Fano lineshapes (see Fig. [2| 

[g n r n /2 + ( g -r/ n )] 2 

ri/4 + (s- Vn ) 2 

where e = A (£i — £2) — 2 7m represents the deviation from 
the n th peak. The asymmetry parameter q n is expressed 
in terms of the two particle scattering matrix 

by the relation 

arg5(27rn/A) = 2g n /(^-l). (7) 

This illustrates the evolution from q n = 00 for free bosons 
(resonance lineshape) to q n —> as c —> 00 (antireso- 
nance). The asymmetry of the lineshape is the first qual- 
itative feature of the crossover regime. The second is the 
finite width T n , for which we give the explicit form below, 
and which vanishes in the two limits. 

The surprising simplicity of our result is a consequence 
of the integrability of the Hamiltonian Eq. Q [9 j. The 
TV-particle scattering it describes is nondiffractive, con- 
sisting of pairwise scattering that either preserves or ex- 
changes the momenta of the scattering particles. A re- 
markable consequence is that the time dependence of the 
A/"-particle propagator describing the amplitude for par- 
ticles at 2/1, • • • , Vn to arrive at xi, . . . xjy after time t can 
be written explicitly for c > as [10] 



(x 1 ,X2,.--,x N \y 1 ,y2,.-.,yN',t)= ^ J "' J A( ^t[ e 
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where Sn denotes the symmetric group of degree TV, and 

K = n {S(K(a) - Ktf)) 'x a <X(3 but y a{a) > y a (/3)}. 

(9) 

To verify Eq. (|| one should first observe that it sat- 
isfies the boundary condition - -^j Gn\ x% =x 3 = 
cQN\xi=xj imposed by the interaction. Next we 



must check that the initial condition £/at(x|v;0) = 
So- Eli $( x i ~ Vcr(j)) is obeyed. This follows from the fact 
that the integral 
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is nonzero only for A a = 1 i.e. when the {xi} are in the 
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FIG. 3. (Top) Diagram representing a term that survives 
integration over Xj, j ^ 1,2. At the top are the positions 
{2/0-1(1)}) an d a t the bottom are {^o- 2 (i)}- The thick red lines 
correspond to the two coordinates x\ and X2 that are not inte- 
grated over in the two-body density matrix Eq. (|3j). (Bottom) 
Usual contribution to the noninteracting HBT effect. 



same order as the {y%}. This in turn is a consequence 
of the following Golden Rule that we will use repeatedly 
for integrals of this type [11: a particle moving to the 
left (right) must be overtaken by another particle moving 
to the left (right). In the present case the Golden Rule 
restricts us to A a = 1, from which the product of 5- 
functions follows. 

To understand the origin of the Golden Rule, consider 



the integral over k a ^ in Eq. (10). The result can be 
viewed as the Fourier transform of the product of factors 
in A a involving fc CT (j), evaluated at Xj — y a (j)i which is 
thus the convolution of the Fourier transforms of these 
factors. Because S(k) is holomorphic in the lower half 
plane for c > 0, its Fourier transform is supported in 
[0, 00). Thus for Xj —y a (j) > we must have at least one 
factor S(k a ^ — k cr ^) where k a ^ appears first (particle is 
overtaken moving to the right). Likewise for xj —y a (j) < 
we must have at least one factor S(k a ^\ — k a ^) where 
fco-y) appears second (particle is overtaken moving to the 
left). 

The time evolution of our array can be found by 
convolving the propagator with the Gaussian initial 
wavepackets 



^iv(x;t) = J g N (x,y-t)Y[(pj(yi)dy 



(11) 



The utility of this expression would seem to be hampered 
by the momentum integrals and the sum over permuta- 
tions in Eq. However, the former may be evaluated 
in the stationary phase approximation at long times 
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where again we have used the variables £j = Xj/t, and 
the stationary phase integral assumes that these are order 
one in the long time limit. In the above A' a denotes 

A ° = II i s ^ a ^) : Xa < x p but y<<*) > y°(p)}' 

(13) 

To evaluate the probability distribution we require the 
'forward and back' propagator 
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In this expression the scattering phases have the explicit form 

A' ai A'; 2 ='[[{S(t a - ^) : ffi(a) > in(J3) but a 2 (a) < <r 2 (/3)} , 
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Unlike the individual A' a , we see that the form of the 
product does not depend upon the ordering of the {xj}. 
Since we need to integrate over all but two of the {xj} 
to find the density correlation function (see Eq. (|3|), this 
fact is tremendously useful, as it tells us that the integrals 



have the same form as Eq. ( 10 ), and allows us to apply the 
Golden Rule. The only non-trivial terms (i.e. without 
<7i(a) = (12(a) for all a) are of the form illustrated in 



Fig. [3] (top). x\ and X2 are exempted from the Golden 
Rule and correspond to the only particles not overtaken. 
Thus we must have 01 (1) = 02(2) and 02(1) = 

Despite this simplification there would still seem to be 



a great many terms to sum in Eq. ( 14 ). We will now show 



that the remaining terms can be grouped according to the 
order of their contribution in the parameter e _2cA , with 
the lower powers amenable to explicit evaluation. Since 
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FIG. 4. Simplest 'new' contribution, smaller by e~ 2cA . 



the parameter cA is the same as the usual Lieb-Liniger 
parameter 7 = c/n, with the density n = A -1 , the use 
of e~ 2cA as a small parameter is not too restrictive. 

Let us first consider the terms that give rise to the 
usual HBT effect in the case of noninteracting particles 
(Fig. [3j bottom). Each of the x a with (J\(a) = 02 (a) 
lying between <ti(1) = 02(2) and 02(1) = 0"i(2) brings 
a factor S(& ~ £a)S(£ a ~ 6) if ^i(l) < <n(2) and 
S(ti ~ Za)S(£a ~ 6) if *i(2) < ^i(l). After integrating 
over the {x a : a ^ 1,2} and convolving with the Gaus- 
sian wavepackets Eq. ( [TT) ), we can sum all such contribu- 
tions in a geometric series to give 
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which generalizes Eq. (|3| to the interacting case. In 
Eq. (16) we have defined the function 



C(£i,6) = 1-2^5(6 -£i -ic) 
^ c -^) 2 erfc(£ [c - z&]) + e^ c+i ^ 2 erfc(£ [c + ^ 2 ]) 



where erfc(x) is the complementary error 

~ l dt. Eq. (16) is shown in Fig. 



erfc(x) = e~ 



function 
2 



In the limit that C(Ci? C2) is close to unity, Eq. (16) can 
be interpreted as a series of Fano lineshapes Eq. (5]) with 
r] n = — Im£, T n = 2(1 — Re£) > 0, and q n as given in 
Eq. ([7]). Fig. [2] (bottom) illustrates the evolution of q n 
between successive peaks from smaller (close to antires- 
onance) to larger values. 

The physical origin of the asymmetry q n lies in the 
scattering phase of particles 1 and 2 with each other, 
while the width T n arises from the collisions of these par- 
ticles with those that they pass, whose momentum has 
a Gaussian distribution and gives rise to a distribution 
of scattering phases. T n vanishes in the limits c — >• 
and c — > 00, but also when t — > 0. In the last case this 
is a consequence of the typical momenta of the particles 
becoming large (except for particles 1 and 2 whose mo- 
menta are fixed by x\ and x^) and the scattering phase 
for their collisions approaching zero (see Eq. Q). 



Let us now show that the remaining contributions are 
small in the parameter e~ 2cA . Consider the first 'non 
HBT' diagram shown in Fig. [4j Evaluating this diagram 
gives the contribution 

!?^ e -2cA e 4iA(£i-£ 2 ) e -£ 2 (e 1 +e 2 ) 

xexp [^ 2 ((c-^) 2 + (c + ^) 2 )] 
x 5(6 - - £1 - ic) 2 S(& - 6 - 2ic) (17) 

where, as always, we ignore the overlap e ~ A2 / A£2 between 
neighboring sites. The two exponential factors e~ cA arise 
from the pole in the upper half plane of xs coming from 
the 5(^3 — £1) factor, and from the pole in the lower half 
plane of £4 coming from the 5(^2 — £4) factor. In the same 
way, one can show that the power of e~ cA appearing in a 
contribution is at least twice the total number of moves 
to the right (or to the left). 

In conclusion, we have shown that the HBT effect of 
interacting particles in one dimension has a number of 
interesting features that distinguish it from the nonin- 
teracting problem, most notably an asymmetry and fi- 
nite width in the peaks of the density-density correlation 
function of the expanded system. The calculation hinges 
upon the integrability of the ID Bose gas, and indeed 
appears to depend essentially upon c > for the form of 
the propagator Eq. ^ to be valid. 
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